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6 ■ ABSTRACT 

We present ultra-high resolution (0.32 km s _1 ) spectra obtained with the 3.9m Anglo- 
Australian Telescope (AAT) and Ultra-High- Resolution Facility (UHRF), of interstellar 
Na I Di, Na I D2, Ca II K, K I and CH absorption toward two high galactic latitude 
^ stars HD 141569 and HD 157841. We have compared our data with 21-cm observa- 

tions obtained from the Leiden/Dwingeloo H I survey. We derive the velocity structure, 
column densities of the clouds represented by the various components and identify the 
clouds with ISM structures seen in the region at other wavelengths. We further derive 
abundances, linear depletions and H 2 fractional abundances for these clouds, wherever 
possible. 

Both stars are located in regions of IRAS 100/mi emission associated with high galac- 
tic latitude molecular clouds (HLCs) : HD 141569 lies, in projection, close to MBM 37 
and the Lynds dark cloud L 134N while HD 157841 is in the vicinity of the MBM 151. 

Toward HD 141569, we detect two components in our UHRF spectra : a weak, broad 
b = 4.5 kms" 1 component at - 15 kms -1 , seen only in Ca II K absorption and another 
component at kms -1 , seen in Na I Di, Na I D 2 , Ca II K, K I and CH absorption. 
The cloud represented by the - 15 kms -1 component, is warm and may be located in 
a region close to the star. The cloud represented by the kms -1 component has a 
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Ca linear depletion 5(Ca) = 1.4 x 10 -4 and shows evidence for the presence of dust, 
consistent with strong 100//m emission seen in this region. The H2 fractional abundance 
/(H 2 ), derived for this cloud is 0.4, which is typically what is observed toward HLCs. 
We conclude that this kms -1 cloud is associated with MBM 37 and L 134N based 
on the presence of dust and molecular gas (CH) and good velocity agreement with CO 
emission from these two clouds. This places HD 141569 beyond MBM 37 and L 134N, 
which are estimated to be at ~ 110 pc. 

In the case of the HD 157841 sightline, a total of 6 components are seen on our 
UHRF spectra in Na I Di, Na I D 2 , Ca II K, K I and CH absorption. 2 of these 6 
components are seen only in a single species. The cloud represented by the components 
at 1.85 kms -1 has a Ca linear depletion S(Ca) = 2.8 x 10~ 4 , indicating the presence of 
dust. The /(H 2 ) derived for this cloud is 0.45 and there is good velocity agreement with 
CO emission from MBM 151. To the best of our knowledge, this 1.85 kms -1 component 
towards HD 157841 is the first one found to have relative line widths that are consistent 
with pure thermal broadening only. We associate the 1.85 kms -1 cloud seen in our 
UHRF spectra with MBM 151 and conclude that HD 157841 must lie beyond ~ 200 pc, 
the estimated distance to MBM 151. 

Subject headings: atomic data — ISM:abundances — ISM:kinematics and dynamics - 
stars:individual (HD 141569 and HD 157841) 
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1. Introduction 

In this paper, we report the results of an ultra- 
high-resolution (0.32 kms -1 ) absorption line study 
of two stars, HD 141569 and HD 157841. Both stars 
were chosen as background sources for detailed spec- 
troscopic studies of high galactic latitude molecular 
clouds (HLCs), as well as gas in the nearby halo. 

HLCs show structure, with dense clumps as small 
as 0.03 pc, which appear to be virially unstable 
(Pound et al. 1990). Some of these clouds appear 
to be embedded in the local hot interstellar medium 
as suggested by the probable detection of X-ray shad- 
owing from two MBM clouds (Burrows and Menden- 
hall, 1991). If this is indeed the case, it is not clear 
what physical process or combination of processes is 
causing such enhanced molecular abundances. This 
has led some authors (Blitz, 1990) to conclude that 
HLCs may represent a completely distinct population 
of clouds. More recent work by Gir et al. (1994) in- 
dicates that HLCs are almost always located along 
filamentary or looplikc HI structures and appear to 
have condensed from atomic gas in situ rather than 
having been entrained in the HI. 

These clouds are an important factor both in the 
composition and the energy balance of the nearby 
interstellar medium. They have molecular masses 
about 50 M Q and the molecular gas contribution of 
these high latitude clouds in the solar neighborhood 
up to 100 pc, is about 5000 M (Magnani 1993). The 
relation and exchange of energy between the molec- 
ular gas in these clouds and the ambient highly ion- 
ized, moderately ionized and atomic gas is unclear 
although it has been suggested by Blitz (1990) that 
compression of these clouds could lead to a phase 
transition of atomic hydrogen into molecular hydro- 
gen. 

Absorption line studies toward stars whose lines of 
sight intercept regions close to HLCs, probe the dif- 
fuse component of the ISM associated with these high 
latitude clouds. Typical values of the kinetic temper- 
atures for the diffuse component of the ISM range 
between 45 to 130 K (Savage et al. 1977). To resolve 
absorption lines arising from clouds with cloud kinetic 
temperatures of = 50 K, a velocity resolution of 
0.30 kms -1 corresponding to a resolving power R = 
10 6 is necessary. The Ultra High Resolution Facility 
(UHRF) at the Anglo Australian Telescope (AAT) is 
currently the world's highest-resolution optical astro- 
nomical spectrograph (Diego et al., 1995, Barlow et 



al., 1995) and is ideal for studies of such cool, inter- 
stellar clouds. 

We have begun a study of the spatial distribution 
and properties of atomic and molecular interstellar 
species towards HLCs at ultra high resolution (veloc- 
ity resolution = 0.32 kms -1 ), and have compiled a 
list of background stars with reliable photometric and 
spectroscopic data and lying in the line of sight to 
HLCs. We have already observed a small fraction of 
this sample, utilizing the UHRF (Blades et al. 1997), 
and in the following sections, we present UHRF ob- 
servations toward HD 141569 and HD 157841. 

1.1. Overview of the HD 141569 sightline 

HD 141569 has been variously classified as an iso- 
lated Herbig Ae/Be star (The et al. 1994) and a 
B9.5V star by Penprase (1992) and an AOVe star by 
Dunkin et al. (1997). It is an emission-line star asso- 
ciated with an IRAS source, IRAS 15473-0346, which 
is not exactly coincident with it. The heliocentric ra- 
dial velocity of the star is -6.4 kms -1 (Frisch, 1987). 
In the direction of HD 141569, which is located in 
galactic coordinates at (I, b) = (4.°2, +36. °9), this 
translates into a LSR radial velocity of -20.1 kms -1 
. The star exhibits excess far-IR emission (Oudmai- 
jer et al. 1992) believed to arise due to emission 
from dust grains present in a circumstellar disk, and 
Sylvester et al. (1996) include it in their list of Vega- 
like systems. 

A high galactic latitude dark cloud complex at (I, 
b) ~ (4°, 36°) lies close to the star in projection. 
This cloud complex includes the Lynds dark clouds 
L 134, L 183 (which is also often referred to as L 134N; 
henceforth we will refer to this cloud as L 134N) and 
L 1780. The dense molecular CO core MBM 36 is 
located within L 134, while the CO core in L 134N 
is referred to as MBM 37 (Magnani et al. 1985). Of 
these two dark clouds, L 134N, centered at (I, b) = 
(6°.0, 36°. 8) and 2° north of L 134, lies closest to 
HD 141569. HD 141569 lies well within the 100/xm 
emission contours of the L 134N dark cloud (Laureijs 
et al. 1991). 

Penprase (1992) estimated a distance of 190±110 
pc for HD 141569 based on spectroscopic and UBV 
and Stromgren photometric data. There are vari- 
ous distance estimates to the L134, L134N and L1780 
complex: 100±50 pc, based on optical extinction and 
surface brightness observations by Mattila (1979), 
110±10 pc based on photometric data by Franco 
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(1989) and 160 pc based on reddenings of a large 
number of stars in this direction by Snell (1981). The 
parallax of HD 141569 from the Hipparcos Catalogue 
is 10.10±0.83 mas, corresponding to a distance of 
99±8 pc. This implies that the distance to the star is 
very similar to that of the L 134, L 134N and L 1780 
complex. Penprase (1992) estimated an E(B-V) of 
0.12 ± 0.09 for HD 141569. The Hipparcos Cata- 
logue lists Tycho magnitudes which transform to give 
(B-V) = 0.078±0.007, which for an A0V star corre- 
sponds to an E(B-V) of 0.10±0.03 (the uncertainty of 
0.03 is due to the uncertainty in the (B-V) calibra- 
tion for dwarf stars). The star falls just outside the 
outermost CO contours which trace MBM 37, the CO 
core in L 134N (Caillault et al. 1995). Based on avail- 
able information, although HD 141569 appears to lie 
approximately at the distance of L134N and clearly 
does not lie within the dark cloud itself, it is not cer- 
tain whether it lies in front of, or behind this dark 
cloud. This point is discussed further in § 3.1.3. and 
our UHRF data indicate that HD 141569 lies beyond 
L 134N and MBM 37. 

1.2. Overview of the HD 157841 sightline 

HD 157841, a B9 star (SIMBAD database) at (I, 
b) = (16°. 8, 15°. 7), lies in projection near the HLC, 
MBM 151, which is located at (I, b) = (21°. 5, 20°.9). 
We have searched the SIMBAD database and related 
literature and have been unable to find any radial ve- 
locity measurement for this star. 

Based on CO velocity agreement, Penprase (1993) 
suggested that MBM 57 located at (I, b) = (5°.l, 
30°.8) is related to MBM 151, which is - 19° away. 
Photometric data of stars in the region around MBM 57 
(Franco, 1989), show an increase in extinction beyond 
100 pc. Penprase (1992) also observed this increase 
in the reddening and estimated a distance of 150 to 
210 pc for the near edge of MBM 57. Assuming that 
MBM 57 and MBM 151 are related, the distance to 
MBM 151 is - 200 pc. The CO (J= l->0) LSR ve- 
locity of MBM 151 is -0.8 kms -1 (Magnani et al., 
1985). The region near HD 157841 has not been 
mapped in CO and it is not known if there is any 
CO emission present at the position of the star or in 
its vicinity. The Hipparcos Catalogue lists a paral- 
lax for HD 157841 of 5.70±1.16 mas corresponding to 
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a distance of 175± 30 pc. It also lists Tycho magni- 
tudes which transform to yield (B-V) = 0.213±0.008. 
For a B9V star, this corresponds to an E(B-V) of 



0.28±0.03. In § 3.2.5., we identify the various compo- 
nents seen in the UHRF spectra toward HD 157841. 

2. The data and its reduction 
2.1. The UHRF data 

The spectra of HD 141569 and HD 157841 were 
obtained during April 25 to 28, 1994 and July 14 
and 15, 1995, using the UHRF at the Coude focus 
of the 3.9-meter Anglo-Australian Telescope. For de- 
tails about the UHRF refer to Diego et al. (1995). 
The configuration of the spectrograph was set to give 
close to the highest achievable resolution which is R 
= 938,000 ± 15,000. We used the 35-slice image sheer 
(Diego, 1993) at the entrance aperture. We masked 
the image sheer down to give 18 slices on the detector 
which allowed for the simultaneous detection of the 
inter-order background on either side of the stellar 
spectrum. 

The detector was a Tektronics CCD (1024 x 1024) 
with a pixel size of 24 x 24 /im. The width of each 
individual slice produced by the image sheer is about 
30 pixels, and the 18 slices filled over half the detector 
format while the order background filled the remain- 
der. The entire CCD was read out with eight times 
binning in the direction perpendicular to the disper- 
sion. 

The instrumental set-up gave a wavelength cover- 
age of 2.5A, requiring separate exposures for each of 
the Na I D lines. Individual exposures were normally 
1800 seconds, and the total integration times were 30 
minutes each for the Di and D2 spectra, 60 minutes 
for Ca 11 K and 180 minutes each for the K 1 and 
CH spectra. Comparison exposures using a Th-Ar 
lamp were used to establish the wavelength scale and 
showed that the plate scale is linear with wavelength 
over the 2.5 A coverage (Barlow et al. 1995). Mea- 
surements of the 6328. 160A stabilized Hc-Ne laser 
line taken during the observing run confirmed that 
we achieved a resolution of 0.32 kms^ 1 (FWHM). 

Standard IRAF reduction packages were used to 
combine individual integrations after cosmic ray events 
and the inter-order background had been removed. A 
scattered light correction of 3.5% was applied to all 
the data, derived from observations of a highly satu- 
rated D line in a different star (HD 24263) recorded 
during the same period. The scattered-light fraction 
for the Ca 11, CH and K I lines is not the same as for 
the Na 1 Di and D 2 lines. In particular for the Ca 11 
K line, inspite of background subtraction (by taking 
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averages of the inter-order regions above and below 
the stellar spectrum), there may still be a residual of 
as much as 10% (Crawford et al. 1994). We do not 
have observations of saturated Ca II K lines as in the 
case of the Na I D lines and therefore do not have ex- 
act values for the scattered light fraction in the Ca II 
line. However, we have compared our total equivalent 
widths with the values published by Penprase (1993), 
assuming that the Penprase data has no scattered 
light contribution. The equivalent widths obtained 
for Na I and CH for both HD 141569 and HD 157841 
are consistent with zero contribution from scattered 
light. For the case of the Ca II line, we estimate the 
contribution of scattered light to the measured equiv- 
alent widths to be - 15% for HD 141569 and - 8% 
for HD 157841. For the case of the K I lines, we are 
not able to make such scattered light estimates since 
there are no published K I equivalent widths for either 
star. We have not subtracted this correction for the 
Ca II lines and estimate that the uncertainty due to 
scattered light has a small effect on the derived col- 
umn densities (< 10%) and has almost no effect on 
the measured linewidths, which is our main concern 
in this paper. 

The spectra were normalized by fits of straight 
lines to sections without telluric or interstellar ab- 
sorption. Telluric absorptions around the D line re- 
gion were removed by dividing by a UHRF telluric 
template spectrum of a Eri, which was rescaled so 
that the depth of individual telluric features matched 
those of unblended telluric lines in the spectra. Rest 
wavelengths in air of 5895.9234 A for Di, 5889.9512 A 
for D 2 , 3933.663 A for Ca n K, 7698.974 A for K I 
and 4300.3133 A for CH, were adopted from Mor- 
ton (1991). The signal-to-noise ratios derived from 
the r.m.s. fluctuation in the stellar continuum for 
HD 141569 arc 20 for D 1; 28 for D 2 and 10 for 
the Ca II K, CH and K I spectra. In the case of 
HD 157841, the S/N values are 35 for Di, 16 for D 2 , 
10 for K I, 20 for Ca n K and 42 for CH. 

2.2. Profile fitting 

We used the method of profile fitting developed 
by Welty et al. (1994) to derive the properties of 
the individual interstellar clouds. The essential as- 
sumptions are that each cloud is homogeneous and 
can be described by a column density N, radial ve- 
locity v, and a Gaussian line of sight velocity dis- 
persion characterized by the line width parameter 
b = (2KT k /m + 2w t 2 ) 1 / 2 , where v t is the r.m.s (Gaus- 



sian) turbulent velocity along the line of sight. The 
intrinsic line profiles were calculated assuming a min- 
imum number of velocity components and were then 
convolved with the instrumental response. The three 
parameters (N, v, b) for each cloud were then deter- 
mined by an iterative, nonlinear, least-square fit to 
the observed profile. The parameters were fine-tuned 
until a satisfactory fit was achieved. The main crite- 
rion of a satisfactory fit is that the r.m.s. deviation 
(tr) from the observed fit is comparable to the r.m.s. 
fluctuation in the stellar continuum. The uncertain- 
ties were determined by varying one parameter while 
keeping others unchanged until the residual between 
the fit and observed profile goes beyond 3a of the stel- 
lar continuum. 

For both HD 141569 and HD 157841, we started 
with the K I profile. To provide a good fit, one compo- 
nent was required for HD 141569 and two components 
for HD 157841. We started with the components de- 
termined from the K I fitting and added more compo- 
nents to generate good fits for the other lines. For the 
case of the Ca II K profile, the prominent features co- 
incide with the components seen in K I but we needed 
to include additional components for both HD 141569 
and HD 157841 to generate good fits. In the case of 
HD 157841 for example, only one of the K I compo- 
nents coincides with the Ca n K components and we 
needed a minimum of 5 more components to achieve 
a good fit. For the Na Di line, the central part is 
saturated for both stars but in the case of HD 157841 
for example, the four components derived from the 
Ca II K fitting, were sufficient to obtain a satisfactory 
fit. The theoretical profile was calculated using the 
same parameters for the Na D 2 line and a satisfactory 
fit was obtained, confirming the choice of v, N and b 
parameters used in the Na Di fit. 

The CH R 2 (1) line is split due to A-doubling in 
the 2 n ground state of the molecule, with the two 
components being separated by 1.43 kms -1 (Black 
& van Dishoeck, 1988). The two A doublet compo- 
nents have equal strength which results in a symmet- 
ric absorption-line profile in the absence of additional 
unresolved velocity structure. The CH line is weak 
toward both stars and the S/N ratios of the spectra 
are probably not enough to clearly see this splitting. 
The slight asymmetry of the profiles suggests the pos- 
sible presence of unresolved cloud structure. We have 
however used only one component to fit the CH profile 
for each of the stars. 
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Penprase (1993) obtained CH spectra for HD 141569 
and HD 157841 at a much lower resolution of A/ S A ~ 
75,000. Although the values of the velocity centroids 
for the CH profiles were not published in this paper, 
a comparison of our UHRF spectra with Penprase's 
spectra shows consistency in the CH velocities. 

Figures 1 and 2 show the normalized profiles plot- 
ted against LSR velocity in the various absorption 
lines for HD 141569 and HD 157841 respectively The 
observed profiles are shown as histograms while the 
best-fit model absorption profiles are shown as contin- 
uous lines. The parameters resulting from the profile 
fitting are listed in Tables 1 and 2 for HD 141569 and 
HD 157841 respectively. 

2.3. H I 21-cm data 

H I 21-cm data from the Galactic survey made 
with the Dwingeloo 25 m telescope (Hartmann 1994, 
Hartmann and Burton 1997) for the regions around 
HD 141569 and HD 157841 in the velocity range -100 
kms -1 to +100 kms^ 1 are presented in Figures 3 and 
4 respectively. The H I observations were made us- 
ing a beam width of 36'. The mean sensitivity of 
the data is about 0.07K and the velocity resolution 
is 1.03 kms -1 . Details of the observations and data 
reduction procedures can be obtained from the two 
references mentioned above. 

The top portions of Figs. 3 and 4 show the posi- 
tions of the HD 141569 and HD 157841 respectively, 
in galactic coordinates, with respect to the positions 
of the H I observations. The lower portions of Figs. 3 
and 4 show H I line profiles at the 9 positions spaced 
0.°5 apart and centered approximately around the 
two stars. We have fitted Gaussian profiles to the H i 
data and estimated the velocities and column densi- 
ties for the individual line components for the regions 
around HD 141569 and HD 157841 which are listed 
in Tables 3 and 4 respectively. Column densities have 
been estimated using equation 8-57 from Mihalas and 
Binncy (1981), assuming that the material is optically 
thin and has a spin temperature of 100 K. As can be 
seen from Figs. 3 and 4, there are large variations in 
line strength from profile to profile. The UHRF ab- 
sorption line profiles are compared to the H I profile 
in Figs. 1 and 2. The top five panels in both figures, 
show the Na Di, Na D 2 , Ca II K, K I and CH absorp- 
tion profiles respectively while H I profile closest to 
the star is shown in the lowest panel. 



3. Discussion 

3.1. The HD 141569 sight line 

As mentioned in § 1.1, the main ISM structures 
in this direction are L 134N and MBM 37. L 134N 
is a prototypical cold, dark cloud with no internal 
energy source and has been extensively mapped at 
millimctric wavelengths. The kinetic temperature in 
the dense region of L 134N has been derived to be 
T k - 12 K by Snell (1981). Swade (1989) has made 
the most detailed mapping of this cloud in the molec- 
ular emission lines of C 18 0, CS, H 13 CO+, SO, NH 3 
and C3H2. The beam sizes used for these observa- 
tions range from 46 to 99" with a grid spacing of ~ 
V and velocity resolution of 0.16 kms -1 . The emis- 
sion lines in all the molecular spectra are relatively 
narrow (FWHM - 0.5 to 2 kms -1 ) and are centered 
around v L sn ranging between 2.16 to 2.65 kms -1 with 
a typical uncertainty of ~ 0.35 kms -1 . These val- 
ues are consistent with the H2CO velocities of ~ 2.4± 
0.7 kms -1 at 2 mm and ~ 2.8±0.9 kms -1 at 2 cm 
obtained by Snell (1981) and the [C 1] and CO veloc- 
ities of ~ 2.3 kms -1 obtained by Stark et al. (1996). 
MBM 37, the molecular core associated with L 134N 
has a similar velocity of ~ 2.3 kms -1 derived from 
the CO(l— >0) observations by Magnani et al. (1985), 
implying close association with L 134N. 

The dark clouds L 134 and L 134N appear to be 
connected to several clouds and have been described 
as resembling "beads on a string" by Magnani et al. 
(1996). These authors report that projected on these 
clouds is an H 1 filament with significantly higher LSR 
velocity. We checked the 21-cm survey of Heiles & 
Habing (1974) and there is a large-scale H 1 feature 
extending between I ~ 10 to ~ 340° and b ~ 36 to 38° 
and mainly at u LSR ~ -3 to +3 kms -1 but extending 
between ~ -10 to +10 kms -1 . There is no clear evi- 
dence for any H 1 at significantly higher LSR velocity 
on the Heiles & Habing plots. Our H 1 profiles at 
position 5 do show a broad component at 6.9 kms -1 
but we are not certain this is the H 1 feature referred 
to by Magnani et al. The column density for this H I 
feature is N(H 1) = 3.53 x 10 20 ± 1.13 x 10 20 cm~ 2 
and it has has no counterpart in the UHRF spectra. 
We conclude that this H 1 component at 6.9 kms -1 
must therefore lie beyond HD 141569. 
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3.1.1. The UHRF spectra 

We base our interpretation of the components seen 
in the UHRF spectra, within the context of the Mc- 
Kee & Ostrikcr (1977) three-phase model. Accord- 
ing to this model, the neutral ISM (the cold neutral 
medium, CNM) is contained in cold (T ~ 80 K), dense 
(n ~ 10 to 1000 cm~ 3 ) clouds which contain most of 
the mass of the ISM but have small filling factors. 
These cold, neutral clouds are surrounded by both a 
warm neutral medium (WNM) and a warm ionized 
medium (WIM). All these phases are embedded in a 
hot (T ~ 10 5 K), highly ionized medium. 

Toward HD 141569, the UHRF spectra show two 
components: one at w LSR = -15.1 kms -1 which is seen 
only Ca n K absorption and another feature at <~ 
kms -1 that shows Na I Di, D 2 , Ca n K, K I and CH 
absorption. These components and their associations 
are discussed in detail in the following sections. 

Despite the ultra-high resolution of 0.32 kms -1 of 
our observations, the b values derived for the 5 compo- 
nents (listed in Table 1), is greater than 1.45 kms -1 
. Welty et al. (1994) made a high resolution survey 
of interstellar Na I Di lines towards 38 stars. Their 
spectra had a resolution of 0.5 kms -1 and an average 
S/N ~ 100. They detected a total of 276 clouds with 
a median b value of 0.73 kms -1 and a median logiV 
value of 1 1 .09 cm -2 . Hence there is the possiblity that 
the b values reported in our paper are overestimates 
of the true b values of individual clouds, caused by 
our inability to resolve the component structure due 
to the relatively low S/N of our spectra. However, 
on the other hand, our earlier UHRF observations on 
the HD 28497 sightline (Blades et al. 1997) have S/N 
comparable to the S/N of the spectra presented here. 
In the case of HD 28497, four of the 15 components 
detected in Na, I had b values ranging between 0.31 
to 0.40 kms -1 while we do not detect any component 
toward HD 141569 (or HD 157841) with b values < 
1.30 kms" 1 . The minimum detectable column den- 
sity for Na I in the spectra presented in this paper 
is ~ 2 x 10 9 cm~ 2 which is about a factor 100 less 
than the column densities of the typical interstellar 
clouds studied by Welty et al. (1994). One possib- 
lity is that the large b value components seen toward 
HD 141569 (and HD 157841) are a superposition of 
several low column density (< 2 x 10 9 cm~ 2 ) clouds, 
with overlapping velocities. These are not the typical 
interstellar clouds described by Welty et al. (1994) 
and it is unlikely that all the components seen in our 



spectra are due to such overlapping low column den- 
sity clouds. For this reason, we conclude that the b 
values we estimate in this paper probably reflect the 
true b values of the individual clouds. 

3.1.2. The -15 kms^ 1 component 

Seen in absorption only in Ca II K, this weak com- 
ponent is relatively broad with b = 4.5 kms -1 . No 
H I counterpart to this component is seen at position 
5. The H I beam (36') is much larger than the ab- 
sorption beam (which is infinitcsimally small) and it 
could be that the H I emission associated with this 
weak component is lost in the general background 
emission. The Ca II column density derived for this 
component is 8.11 x 10 11 cm~ 2 . Unfortunately, in- 
terstellar Ca II column densities are very poorly cor- 
related with quantities like N(H i), N(Na i), N(K i) 
and E(B-V) (Hobbs, 1974) and it is not a useful ex- 
ercise to derive N(H i) using the observed N(Ca n) 
for this component. We can put an upper limit to 
the N(H i) of this component, corresponding to the 
sensitivity of our H I data which is ~ 10 18 cm~ 2 . As- 
suming that Ca II is the dominant ionization stage of 
this element in the gas, this implies the linear deple- 
tion 5(Ca) - [N(Ca)/N(H)]/[N(Ca)/(N(H)] Q > 0.12, 
implying relatively low depletion for this component. 

In the absence of turbulence, b = 4.5 kms -1 for 
this component, corresponds to a kinetic temperature 
< 47,900 K and such a warm component may be 
located in a region close to the star. HD 141569 has 
a w LSR radial velocity of -20.1 kms -1 and the asso- 
ciation of the -15.1 kms -1 component with a region 
close to the star is not clear. 

3.1.3. The kms^ 1 component 

The Na I Di, D 2 velocities for this component 
match the Ca II K, K I and CH component veloci- 
ties to better than 0.32 kms -1 . In the absence of 
turbulent motions, the kinetic temperature of an ab- 
sorbing component for Na I is given by = 1393 b 2 , 
for Ca II by T k = 2422 6 2 , for K I by T k = 2368 b 2 
and for CH by T k = 787 b 2 . If all these absorption 
components arise from coextensive gas where turbu- 
lent motions are absent, the b values for the various 
species should differ systematically : 6(Ca) should 
be 0.756(Na), 6(K) should be 0.77&(Na) and 6(CH) 
should be 1.336(Na). From Table 1, it is clear that 
this is not the case since the b values for Ca II and 
K i exceed the &(Na) value and the 6(CH) value is 
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> 1.33 x £>(Na). On the other hand, if turbulent mo- 
tions were dominant in the region where these absorp- 
tion components arise, the b values for all the species 
should be exactly the same and this is also not what 
is observed. Therefore the various components must 
arise from different regions of a cloud, which repre- 
sent different phases of the ISM. This is similar to 
the interpretations by Barlow et al. (1995) for clouds 
toward £ Oph and Blades et al. (1997) for clouds 
toward HD 28497. 

In the absence of turbulent motions, (Na) = 
2940 K and T k (Ca) = 8570 K and we interpret the 
Na I component to arise from gas associated with the 
cold neutral medium (CNM), enveloped by the warm 
gas delineated by the Ca component. The K I com- 
ponent has a T^. = 12,310 K. It is possible that the 
K I absorption feature has velocity structure but the 
S/N ratio of 10 for the spectra, has led us to estimate 
a large b(K) for this weak absorption feature. This 
issue can be resolved by higher S /N UHRF spectra in 
K I. If these future higher S/N UHRF spectra, still 
show no velocity structure for this weak K I feature, 
it would imply that a warm K I component envelopes 
the cooler Ca n K and Na I components. 

The kinetic temperature obtained for the CH com- 
ponent is Tjj. = 3910 K, assuming no turbulent mo- 
tions. The CH feature is weak and the S/N ratio of 
42 for the CH spectra is probably not sufficient to re- 
solve the two components of the CH molecule and the 
b value that we have obtained may include a blend of 
both components, resulting in a large value. CH 
spectra at higher S /N obtained with the UHRF would 
be required to clarify this issue of additional velocity 
components. If the higher S/N spectra still show no 
velocity structure, then the CH component toward 
HD 141569 must be associated with the warm, outer 
halo surrounding the cold molecular core of MBM 37. 
A similar interpretation was made for the broad b = 
1.8 kms" 1 CH component toward £ Oph by Craw- 
ford et al. (1994). The CO gas associated with both 
L 134N and MBM 37 (refer § 3.1) is at w LSR - 2.4 
kms -1 and association of the kms -1 component 
with MBM 37 would imply that HD 141569 lies be- 
yond ~ 110 pc. 

We are interested in determining abundances and 
depletions for the entire cloud represented by the com- 
ponents at ~ kms -1 seen in the various species. 
Although the various absorption components arise in 
different regions of the cloud, it is valid to use the col- 
umn density values derived from the UHRF spectra 



for the various species for the purposes of studying 
abundances of the ~ kms -1 cloud. Based on ve- 
locity agreement, we identify the H I component at 
position 5 with t> LSR = 0.78 kms -1 to be associated 
with this kms -1 UHRF component. This compo- 
nent has N(H i) = 5.03 x 10 20 ± 5.57 x 10 18 cm" 2 
(Table 4). The linear depletions derived are 5(Na) 
= [N(Na)/N(H)]/[N(Na)/(N(H)] Q = 0.004, 5(Ca) = 
[N(Ca)/N(H)]/[N(Ca)/(N(H)] = 1.4 x 10~ 4 and 
S(K) = [N(K)/N(H)]/[N(K)/(N(H)] = 0.0017 im- 
plying extremely high depletions for this cloud. High 
depletions imply low gas phase abundances and there- 
fore high dust phase abundances since the missing 
elements are believed to be incorporated into dust, 
particularly in the case of calcium (Barlow & Silk, 
1977, Barlow 1978). Evidence for the presence of dust 
in this component from the UHRF data is consistent 
with the strong 100/um emission from L 134N (Lau- 
reijs et al. 1991). 

The H2 column densities in diffuse interstellar 
clouds are well correlated with the CH column densi- 
ties through the relation N(H 2 ) = 2.6 x 10 7 N(CH) 
with an r.m.s. scatter of 9.0 x 10 19 (Somcrvillc & 
Smith, 1989). Using this relation and N(CH) derived 
from the UHRF data, we obtain N(H 2 ) = 1.51 x 
10 20 ± 9.0 x 10 19 cm- 2 . The H 2 fractional abun- 
dance, given by /(H 2 ) = 2N(H 2 )/[N(H 1) + 2N(H 2 )], 
for this kms -1 cloud is therefore ~ 0.4±0.2. This 
/(H 2 ) value which is typically what is observed for the 
diffuse component associated with HLCs (Penprase, 
1993). Both the evidence for the presence of dust 
and molecular gas and the H 2 fractional abundance 
value for this kms -1 cloud strengthen the case for 
its association with MBM 37 and L 134N. 

3.2. The HD 157841 sight line 

A total of 6 components are seen toward HD 157841 
in Na I, Ca 11 K, K 1 and CH (Fig. 2 and Table 
2) in the UHRF absorption spectra. The absorption 
features in the various species broadly resemble each 
other, though they differ in details. 2 of these 6 com- 
ponents are seen only in a single species. A com- 
parison of the UHRF spectra with the H 1 profile at 
position 4, located closest to HD 157841 (lowest panel 
of Fig. 2) shows that the -8.8 kms -1 and 1.85 kms -1 
UHRF components have H 1 counterparts. There is 
an additional H 1 feature at 14.5 kms -1 which has no 
corresponding absorption feature on the UHRF spec- 
tra and this H 1 emission must be associated with 
gas beyond HD 157841. In the following sections, we 
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discuss the individual components and their nature. 

3.2.1. The -15 fans" 1 component 

This weak, broad b = 4.97 kms" 1 component is 
seen only in Ca n K absorption. Assuming no tur- 
bulent motions, the kinetic temperature derived for 
this component is = 59,900 K and this compo- 
nent is likely to be associated with a region close to 
the star. As mentioned in § 1.2, we have been unsuc- 
cessful in finding any stellar radial velocity determi- 
nations for HD 157841 and therefore we cannot make 
any further comments on the association of this cloud 
with a region close to the star. This Ca II component 
has no associated H I emission which could be due 
to the larger H I beam size causing weak components 
with small angular extent to dissappear in the general 
emission. 

3.2.2. The -8.8 km s^ 1 component 

This component is seen in Na I and Ca n K ab- 
sorption with 6(Na) ^ 0.756(Ca) or 6(Na) = &(Ca) 
indicating that they do not arise from coextensive 
gas. Assuming no turbulent motions, (Na) = 
2370 K and T k (Ca) = 6595 K. We interpret the 
observed warm Ca n component to envelope the 
colder Na I cloud core, similar to what is observed 
in case of C Oph (Barlow et al. 1995), HD 28497 
(Blades et al. 1997) and the kms" 1 component 
toward HD 141569 (§ 3.1.3.). The H I emission as- 
sociated with this component has N(H i) = 4.89 x 
10 19 ± 5.53 x 10 18 cur 2 . This results in 5(Na) = 
[N(Na)/N(H)]/[N(Na)/(N(H)] = 0.0078 and <5(Ca) 
= [N(Ca)/N(H)]/[N(Ca)/(N(H)] Q = 0.002, indicating 
relatively high depletions and implying the presence 
of dust. 

3.2.3. The -1.95 fcms" 1 component 

This component is seen in Ca II K and K I absorp- 
tion. The Ca n feature is relatively narrow with b = 
1.68 kms -1 and assuming no turbulent motions, the 
kinetic temperature T k (Ca) = 6820 K. The K I fea- 
ture is narrow with b = 0.71 km s" 1 which corresponds 
to a Ty. (K) = 1210 K in the absence of turbulence. 
The warm Ca II K component is interpreted as the 
warm envelope surrounding a cooler K I component. 

3.2.4- The 0.66 fcms" 1 / 1.85 kms -1 component 

The 1.85 km s -1 component is a strong feature seen 
in the UHRF spectra in Na I, Ca II and K I absorp- 



tion. As noted in § 2.2, the CH absorption line toward 
HD 157841 has a slight asymmetry and suggests pos- 
sible unresolved cloud structure. We believe that due 
to this asymmetry, this 0.66 kms 1 component seen 
only in CH is associated with the 1.85 kms -1 com- 
ponent seen in Na I, Ca n and K I. 

For the Na I, Ca n and K I absorptions, &(Na) « 
0.756(Ca) (see Table 2) w 0.776(K). The systematic 
difference in the linewidths in these species implies 
they arise from coextensive material where turbulent 
motions are not important. The kinetic tempera- 
tures corresponding to the derived b- values of the 1.85 
kms -1 component are 5630 K for Na I, 5230 K for 
Ca II and 6370 K for K I, in good agreement with each 
other. This component therefore appears to corre- 
spond quite well to the warm neutral medium (WNM) 
phase posited by the McKee and Ostrikcr model of the 
ISM. The CH absorption has &(CH) = 3.19 kms -1 
which is > 1.336(Na) (refer § 3.1.3.). However, as in 
the case of the CH spectra of HD 141569, the low 
S/N of the spectra may have led us to overestimate 
the b- value for this component. 

The velocity of this component closely matches the 
CO LSR velocity of MBM 151, which is - 0.8 kms" 1 
. This points to an association between this compo- 
nent and MBM 151 and places HD 157841 beyond 
MBM 151 which is at - 200 pc (see § 1.2). In the 
absence of turbulence, the kinetic temperature 
(CH) derived for this component is ~ 8000 K and this 
molecular gas may be a warm envelope surrounding 
the cold core of MBM 151. 

There is a strong H I 21-cm emission at 2.49 kms" 1 
which is most likely associated with the UHRF fea- 
ture. The small velocity difference between this H i 
feature and the Na I, Ca n and K I absorptions is not 
significant if one takes into consideration the differ- 
ences in beam sizes for the H I and UHRF observa- 
tions. The column density of this feature is N(H i) = 
1.27 x 10 21 ± 2.15 x 10 19 cm" 2 . Using the Somcrvillc 
& Smith (1989) relation (§ 3.1.3.) and the UHRF 
derived CH column densities, we estimate N(H2) = 
5.02 x 10 20 ± 9 x 10 19 cm" 2 . The H 2 fractional 
abundance, /(H 2 ), for this cloud is 0.45±0.08 which 
is typically what is found in HLCs (Penprase, 1993) 
and which is similar to the value determined earlier 
in § 3.1.3. for the kms" 1 component seen toward 
HD 141569. 

The linear depletions derived for this cloud are 
<S(Na) = [N(Na)/N(H)]/[N(Na)/(N(H)] = 0.023, <5(Ca) 
= [N(Ca)/N(H)]/[N(Ca)/(N(H)] = 2.8 x 10" 4 and 
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5(K) = [N(K)/N(H)]/[N(K)/(N(H)] Q = 0.0044 for 
this component. This cloud is therefore highly de- 
pleted and has high dust phase abundances. Both 
the evidence for the presence of dust and the H2 frac- 
tional abundance values, make the association of the 
0.66 kms -1 CH component and the 1.85 kms -1 Na I, 
Ca II and K 1 components with the same cloud and 
with MBM 151, more secure. 

3.2.5. The 4-65 kms^ 1 component 

This component is seen in both Na 1 and Ca 11 ab- 
sorption. 6(Na) w 0.756(Ca) for this component and 
the two features must arise from coextensive material. 

4. Conclusions 

We have made an ultra-high resolution study of 
the Na I, Ca II K, K I and CH interstellar absorption 
lines toward two stars, HD 141569 and HD 157841. 
These absorption spectra have been compared to 21- 
cm data obtained from the Leiden/Dwingeloo H 1 sur- 
vey. Both stars probe the gas in regions close to high 
galactic latitude molecular clouds : the HD 141569 
sightline intercepts a region close to MBM 37 and 
L 134N while the HD 157841 sightline probes a re- 
gion close to MBM 151. 

The results of our investigation are as follows. 
1) Toward HD 141569, two components are seen in 
our UHRF spectra : one at -15 kms -1 and another 
at kms" 1 . The -15 kms -1 component, seen only 
in Ca II K, is weak, broad and shows low linear deple- 
tion, <5(Ca) > 0.12. In the absence of turbulence, the 
kinetic temperature derived is < 47,900 K and the 
cloud represented by this component may be located 
in a region close to HD 141569, whose stellar radial 
LSR velocity is -20.1 kms -1 . 

The ~ kms -1 component is seen in Na 1, Ca 11 
K, K 1 and CH absorption. The cloud represented by 
these components has a linear Ca depletion of <5(Ca) = 
1.4 x 10~ 4 , implying the presence of dust, consistent 
with strong 100 /im emission from this region, and H2 
fractional abundance /(H 2 ) = 0.4. We conclude that 
this kms -1 cloud is associated with MBM 37 and 
L 134N because of the presence of dust and molecular 
(CH) gas at velocities close to CO velocities derived 
for the two HLCs. The Hipparcos distance of 99±8 pc 
for HD 141569 is comparable to the distance of 100 - 
110 pc that has been estimated for the MBM 37 and 
L 134N clouds. However, our absorption line data 
imply that these clouds must lie a little closer than 



HD 141569. 

2) A total of 6 components are seen in our UHRF 
spectra in Na I, Ca 11 K, K 1 and CH absorption to- 
ward HD 157841. The cloud represented by the 1.85 
kms -1 absorption components has a Ca linear deple- 
tion <5(Ca) = 2.8 x 10~ 4 and the derived H2 fractional 
abundance is /(H2) = 0.45. The relative ^-values de- 
termined for the 1.85 kms -1 Na 1, Ca 11 and K 1 
absorption lines are found to be consistent with ther- 
mal broadening (with no significant turbulent compo- 
nent), corresponding to a kinetic temperature of 7^ 
= 5700±500 K. To our knowledge, this is the first ve- 
locity component found whose line widths can unam- 
biguously be attributed to pure thermal broadening. 
We conclude that the 1.85 kms -1 cloud is associated 
with MBM 151 because of the presence of dust and 
molecular (CH) gas at velocities close to CO veloci- 
ties derived for MBM 151. HD 141569 must therefore 

lie beyond MBM 151, estimated to be at ~ 200 pc 

45 

and its Hipparcos distance of 175±3o pc is consistent 
with this conclusion. 

This research has made use of the SIMB AD database, 
operated by the CDS Strasbourg, France. MSS was 
supported by grant GO-6723 from Space Telescope 
Science Institute. 
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Table 1 

Model Fit Parameters for HD 141569 



Na I Di, D 2 C a II K K I CH 



f LSR a ••• -15.1 c 

iV b ••• 2.83 d 

6 a ••• 4.5 d 

w LSR a 0.07 c -0.11 c 0.23 c 0.09 

iV b 45.4 C 1.68 c 1.16 d 58.0 C 

& a 1.45 c 4.45 c 2.28 d 2.23 c 



a w LSR and b arc in units of km s 
b iV is in units of 10 n cm" 2 
c crror: < 10% 
d error: 10 to 30% 
e error: > 30% 
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Table 2 

Model Fit Parameters for HD 157841 



Na I Di, D 2 C a ii K K i CH 
W 1 ••• -14.97 d 







1 QA d 

J. . -J I 








4 97 d 




^LSR a 


-8.89 c 


-8.7 C 




iV b 


7.81 c 


2.40 c 






1.30 c 


1.65 c 








-2.17 c 


-1.95 c 


iV b 




2.09 c 


2.39 c 


6 a 




1.68 c 


0.71 c 










iV b 
















^LSR a 


*1.85 c ' f 


1.85 c 


1.85 c 


iV b 


605. C 


841 c 


7.53 c 




2.01 c 


1.47 c 


1.64 c 




*4.65 c 


4.65 c 




iV b 


5.78 c 


1.88 c 




6 a 


1.91 c 


1.52 c 





0.66 c 
193.0 C 
3.19 c 



a t> LSR and b arc in units of kms 

b 7V is in units of 10 n cm~ 2 

* component not well resolved in the spectra 

^possible weak counterpart to the Ca II and K I at - 
1.95 kms -1 is hidden in this saturated Na I component 

c error: < 10% 

d error: 10 to 30% 

°error: > 30% 
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Table 3 

Model Fit Parameters to the HI line profiles for HD 141569 



Position no. 


I, b 


velocity (km s : ) 


FWHM (km s^ 1 ) 


N(HI) 10 19 cm" 2 


1 


3.°5, 36. °5 


-3.14±0.53 


15.03±0.74 


1.06±0.06 






0.69±0.01 


5.48±0.03 


9.39±0.07 






13.52±0.31 


14.27±0.31 


1.14±0.02 


2 


3.°5, 37.°0 


-0.26±0.96 


19.58±1.21 


1.09±0.08 






n kq+ci ni 

u . ou znu . u i 


c OQJ-n HQ 
O .iOXU.UO 








15.18±0.86 


16.21±0.93 


0.86±0.09 


3 


3.°5, 37.°5 


0.46±0.01 


4.77±0.03 


9.74±0.05 






4.65±0.53 


24.05±1.07 


1.42±0.03 






20.25±0.61 


9.98±2.78 


0.22±0.07 


4 


4.°0, 36. °5 


-3.08±0.47 


14.09±0.80 


17.30±1.40 






0.81±0.01 


5.43±0.04 


54.36±0.69 






13.19±0.40 


15.45±0.76 


17.84±1.05 


5 


4.°0, 37.°0 


0.78±0.02 


4.75±0.04 


50.3±0.6 






6.92±0.39 


25.63±0.80 


35.3±1.13 


6 


4.°0, 37.°5 


0.54±0.01 


4.99±0.02 


47.98±0.24 






3.80±0.38 


17.69±0.75 


11.84±1.35 






9.66±0.48 


31.52±0.75 


23.20±1.37 


7 


4.5, 36. °5 


-3.14±0.65 


14.52±0.92 


17.56±1.94 






0.92±0.01 


5.35±0.04 


53.42±0.75 






12.80±0.60 


16.45±0.95 


17.70±1.48 


8 


4.°5, 37.°0 


0.73±0.01 


5.20±0.03 


51.37±0.50 






3.65±1.25 


23.33±1.80 


30.62±3.77 






18.13±3.44 


19.84±4.58 


5.25±3.80 


9 


4.°5, 37.°5 


0.60±0.01 


4.86±0.03 


45.07±0.43 






2.88±0.44 


22.19±1.11 


27.24±0.88 






20.10±0.50 


12.17±1.27 


4.71±0.95 
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Table 4 

Model Fit Parameters to the HI line profiles for HD 157841 



Position no. 


J U 

I, 


velocity (km s ) 


r Wxllvl (km s J 


IN (til) 1U cm 


1 


16. °5, 15. °5 


-0.42±0.12 


5.81±0.14 


43.70±2.81 






5.04±0.13 


6.99±0.18 


63.53±3.05 






5.89±0.19 


27.84±0.56 


51.45±0.98 


2 


16. °5, 16. °0 


-0.99±0.05 


3.30±0.15 


10.87±0.98 






3.00±0.06 


8.54±0.14 


81.44±2.04 






5.85±0.36 


29.79±1.03 


44.77±1.53 


3 


16. °5, 16. °5 


-9.75±0.39 


6.65±1.06 


6.32±0.78 






0.81±0.24 


6.79±0.36 


56.60±11.78 






5.03±0.99 


9.14±0.74 


40.40±10.93 






18.78±1.81 


20.02±3.44 


17.33±3.29 


4 


17. °0, 15. °5 


-8.74±0.4 


4.86±0.63 


4.82±0.54 






2.49±0.10 


9.89±0.19 


125.0±1.99 






14.51±1.28 


22.92±2.20 


21.7±2.54 


5 


17. °0, 16. °0 


-8.89±0.18 


4.96±0.47 


4.37±0.50 






0.91±0.10 


6.59±0.15 


62.11±7.18 






4.50±0.45 


8.57±0.32 


44.74±6.94 






9.01±0.81 


36.78±1.35 


32.44±1.71 


6 


17.°0, 16. °5 


-8.91±0.26 


4.79±0.63 


33.18±0.53 






1.28±0.08 


6.59±0.12 


69.41±4.73 






5.51±0.75 


8.84±0.58 


20.62±4.49 






10.42±1.00 


38.42±1.67 


30.98±1.87 


7 


17°5, 15. °5 


-9.27±0.12 


3.85±0.30 


2.90±0.27 






0.22±0.10 


6.02±0.15 


35.14±5.18 






3.66±0.20 


7.88±0.16 


72.82±5.24 






7.99±0.35 


31.54±0.63 


40.64±1.04 


8 
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Fig. 1. — Normalized profiles (histograms) plotted 
against LSR velocity for HD 141569. The top 5 pan- 
els show UHRF absorption line data in the Na Di, 
Na D2, Ca 11 K, K 1 and CH lines respectively. The 
best-fit model absorption lines profiles are plotted as 
continuous thin lines. The lowest panel shows the H I 
profile at position 5 (refer Fig. 3 and § 2.3), the obser- 
vation which is closest to the location of HD 141569. 
Note the presence of the strong absorption feature at 
kms" 1 which coincides with the H I peak. This 
component shows the presence of both molecular as 
well as atomic hydrogen. For details see § 3.1.3. 

Fig. 2. — Normalized profiles (histograms) plotted 
against LSR velocity for HD 157841. The top 5 pan- 
els show UHRF absorption line data in the Na Di, 
Na D2, Ca 11 K, K 1 and CH lines respectively. The 
best-fit model absorption lines profiles are plotted as 
continuous thin lines. The lowest panel shows the H I 
profile at position 4 (refer Fig. 4 and § 2.3), the obser- 
vation which is closest to the location of HD 157841. 
Note the presence of the strong absorption feature in 
the Na I, Di, D 2 , Ca 11 K and K 1 lines at 1.85 kins -1 

Fig. 3. — The top part of this figure shows the po- 
sition of HD 141569 with respect to the positions of 
the beams used for the nine H 1 observations presented 
in § 2.3, in Galactic coordinates. The asterisk indi- 
cates the position of the star which falls just within 
the beam of position 5. The lower part of the Figure 
shows the H 1 profiles plotted against LSR velocity, 
at the 9 positions taken from the Hartmann & Bur- 
ton (1995) survey. Parameters of Gaussian fits to the 
profiles are listed in Table 3. 

Fig. 4. — The top part of this figure shows the posi- 
tion of HD 157841 with respect to the positions of the 
beams used for the nine H 1 observations presented 
in § 2.3, in galactic coordinates. The asterisk indi- 
cates the position of the star which falls just within 
the beam of position 4. The lower part of the Fig- 
ure shows the H 1 profiles in the LSR velocity range 
-100 to +100 kms -1 , at the 9 positions taken from 
the Hartmann & Burton (1995) survey. Parameters 
of Gaussian fits to the profiles are listed in Table 4. 
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